Bistability and hysteresis in tilted sandpiles 
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We show that tilting a model sandpile that has dynamic 
disorder leads to bistability and hysteresis at the angle of 
repose. Also the distribution of local slopes shows an inter- 
esting dependence on the amount of tilt - weakly tilted sand- 
piles retain the quasi-continuous distributions of the steady 
state, while large tilt makes the distribution more discrete, 
with local slopes clustered round particular values. These ob- 
servations are used to explain recent experimental results on 
avalanche shapes; we give a theoretical framework in terms of 
directed percolation. 

05.50.+q, 45.70.-n, 81.05.Rm 

The angle of repose of a sandpile, 6 r , is the maximum 
inclination of the free surface of a stationary pile. It is 
well known jjj that sandpiles exhibit bistable behaviour 
at and around this angle; this corresponds to a range of 
values for the measured angle of repose which varies as a 
function of different configurational histories. It is con- 
ventional to define this range in terms of another angle, 
m , called the angle of maximal stability; this is the mini- 
mum value of the angle of the sandpile at which avalanch- 
ing is inevitable. Clearly, 9 m > 9 r and the range rn - r 
corresponds to a range of angles between the free sand- 
pile surface and the horizontal such that either a stable 
stationary state or avalanching can result depending on 
how the sandpile was produced. 

In this Letter, we show that such bistable behaviour 
is obtained when we 'tilt' a model sandpile that has dy- 
namic disorder Our findings on the correlation be- 
tween avalanche shapes, and the angle of tilt of the under- 
lying sandpile surfaces, match recent experimental results 
H ; additionally, we provide a theoretical explanation for 
our results in terms of directed percolation. 

Our model sandpile is a two-dimensional version of 
an earlier disordered and nonabelian cellular automaton 
model ^j. The 'grains' are square prisms with dimen- 
sions 1 x 1 x a, and these are placed on sites i,j of a 
square lattice with 1 < i,j < L, where L is the system 
size. A grain within column i,j may rest on either its 
square (lxl) face or a rectangular (1 x a) face. We de- 



note these two states pictorially by — or | because they 
contribute respectively Az = a, 1 to the total columnar 
height z(i,j). 

Grains are deposited on the sandpile with a given prob- 
ability of landing in the — or the | orientation. The square 
face down (— ) configuration of grains is considered to be 
more stable and this implies that in general, and certainly 
well away from the surface, grains contribute Az — a to 
the column height. However incoming grains, as well as 
all other grains in the same column, can 'flip' to the other 
orientation with probabilities: 

P{ ► |) = exp(-d/d-) 

P{\ -> -) = exp(-d/d\) (1) 

where d-, d\ are scale heights. This 'flip' embodies the 
elementary excitation involved in the collective dynamics 
of clusters since, typically, clusters reorganise by grain re- 
orientation. The depth dependence reflects the fact that 
surface deposition is more likely to cause cluster reor- 
ganisation near the surface, than deep inside the sand- 
pile. After deposition and possible reorganisation, each 
column has a local slope s(i,j) given by: 

s(i,j) = z(i,j) - -(z(i + + z(i, j + 1)) (2) 

If s(i,j) > s c , where s c is the critical slope, then 
the two uppermost grains fall from column i,j onto its 
neighbours Q or, when i or j = L, exit the system. 
This process could lead to further instabilities and hence 
avalanching. 

Since the critical slope is the threshold for permissi- 
ble slopes, it is a strong determinant of the allowable 
granular configurations in the sandpile, and thus a deter- 
minant of the relative populations of the — and | states. 
It is plausible that s c is inversely proportional to the 'tilt' 
(the angle made by the base of a sandpile with the hor- 
izontal). The greater the tilt the more unstable will be 
the pile to deposition, with a consequent decrease in s c 
(the critical height for column height differences). 

In the following, we fix a, d- and di, and examine the 
effect of varying s c on the mean slope < s > of the sand- 
pile: 
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This slope corresponds to the macroscopic slope taii9 r 
measured in experiments; its variation with s c is shown 
in Fig. |l| as a line of crosses. We note that tan6> r de- 
creases proportionately with s c . Since the free surface 
of the sandpile retains an approximately constant angle 
with the horizontal, tilting it through larger angles (i.e. 
decreasing s c ) would result in lower measured angles of 
repose 8 r , as observed. 

Next, we mimic the effect of suddenly tilting a sandpile 
constructed with a critical slope of s c by reducing this to 
some s' c < s c . A set of stable slopes becomes unstable, 
avalanching occurs, and the sandpile stabilises to a new 
mean slope < s > - except where the angle of tilt is so 
large that spontaneous flow occurs. (Spontaneous flow, 
shown by the full line in Fig. [l] sets in for critical slopes 
s' c < 2a, since such thresholds make even ordered stack- 
ings of flat, i.e. ' — ', grains, unstable). Also, in Fig. [l], 
we show by triangular and circular symbols (correspond- 
ing to triangular and uphill avalanches 0) the nature 
of avalanche 'footprints' generated when a steady-state 
sandpile is tilted through different angles and then per- 
turbed by additional deposition. For example, when a 
sandpile built with s c = 2.05 is tilted so that s' c ~ 1.75, 
avalanches generated by further deposition are, on aver- 
age, triangular in shape (Fig. ||a). Beyond this value of 
s' c , uphill avalanches result (Fig. ^|b). Our results are in 
qualitative agreement with recent experiments || . 

Also, we observe that sandpiles constructed at a given 
tilt (with s c — 1.85, say) have a lower angle of repose 
than those constructed with smaller tilt (corresponding 
to, say s c — 2.05) that are subsequently tilted to the 
same angle; in this case the respective values of the mean 
slope are < s >= 1.21 and 1.39. Since the angle of re- 
pose represents the threshold for stationary behaviour in 
a sandpile, this is a strong indicator of bistable behaviour 
; additionally, hysteresis is manifest in these results since 
they show that the sandpile can either be stationary or 
avalanching, depending on its configurational history, for 
tan6» r e [1.21, 1.39]. 

Hysteresis also manifests itself in the distributions of 
local slopes for tilted piles (Fig. ||). Fig. |^a shows the 
smooth distribution of local slopes obtained when a pile 
is constructed in the steady state with s c = 1.85. Fig. ||b 
shows the much more discontinuous values of local slopes 
obtained when after construction with s c = 2.05, the 
sandpile is tilted to s' c — 1.85. Finally, Fig. |^c shows the 
extremely jagged distributions of local slopes obtained for 
even larger tilt s' c = 1.75. Clearly, different preparation 
histories of our model sandpile in the same macroscopic 
state can lead to extremely different surface conforma- 
tions, as in reality 0. 

The observation that steady-state sandpiles have a 
smooth distribution of local slopes, while tilted sandpiles 
manifest a clustering of values of local slopes is easily ex- 
plained. The upper layers of a steady-state sandpile sub- 



ject to sudden tilt, are 'avalanched' away, exposing bulk 
grains to form the new surface. These tend to be in the 
— orientation, with Az = a; hence slope differences tend 
to be clustered around multiples of a/2. (This is in clear 
contrast to the situation at the surface of steady-state 
sandpiles, where | and — orientations are stochastically 
generated due to deposition and rearrangement). The 
effect of tilting the sandpile through larger angles is to 
bring grains to the surface which were ever deeper in the 
bulk; these tend to be more and more ordered as a func- 
tion of their depth leading to increasingly discrete 
distributions of local slopes. 

Only minor changes in the continuous distribution of 
local slopes occur in the limit of small tilts. An avalanche 
generated by deposition in these circumstances will prop- 
agate only a small distance up the slope (typically after 
traversing only a few sites) until a local slope is encoun- 
tered for which the movement of a grain is insufficient to 
achieve the toppling threshold. The resulting avalanches 
are 'triangular' here, slope instabilities propagate 
downwards most of the time. For sandpiles tilted through 
large angles, on the other hand, only a few discrete val- 
ues of local slopes are permissible (Kig.^Jc) ; the propa- 
gation of grains disturbs the stability of grains uphill of 
the initiation site with a greater probability, causing the 
propagation of 'large' (|] or 'uphill' || avalanches. 

Two possible measures of these distinct avalanche mor- 
phologies are shown in Fig. ||, as a function of s' c . The 
crosses denote the average fraction of the avalanche foot- 
print that is uphill from its point of initiation iq, jo, while 
the dots denote the average distance of the centre of the 
avalanche below the point of initiation. We see that there 
is a sharp transition, as shown by both indicators, be- 
tween triangular and uphill avalanches at s' c ~ 1.8. For 
lower tilt angles (higher s' c ), the avalanche centres are 
well below the point of initiation, and the fraction of 
the avalanche uphill from this is small. For higher tilt 
angles (lower s' c ), the avalanche centres are much closer 
to the initiation point, thus revealing that a substantial 
fraction of the avalanche is uphill H . Similar data, for 
other values of the initial base angle s c , were used to 
construct the 'phase diagram' in Fig. [j] distinguishing 
the two avalanche morphologies. It is easy to see why 
no (discrete) uphill avalanches are observed for sandpiles 
prepared with initially large tilt (s c « 1.9 in Fig. [I]); 
further large tilts simply bring the sandpile to its maxi- 
mal angle of stability with respect to the horizontal, 9 m . 
Denoted by the full line in Fig. [j], this is the angle at 
which a stationary sandpile becomes unstable; sponta- 
neous flow sets in, with continuum avalanches resulting 
from the destabilisation of all sites. 

Differences in the local slope distributions before and 
after tilting, (Fig. 0) have implications for surface rough- 
ness in sandpiles [gUj]] . The quasi-continuous distribution 
of local slopes (Fig. ||a) of steady-state sandpiles result 
in a more continuous, locally 'smoother', surface with 
a large number of disordered grain configurations. On 
the other hand a more discontinuous, locally 'rougher', 
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surface results when there are very few permissible local 
slopes (Fig. ||c) as in the case of the nearly ordered con- 
figurations exposed when sandpiles are strongly tilted. 
However, it appears H that the disordered surfaces of 
weakly tilted sandpiles are globally rougher, despite be- 
ing locally smoother in the sense of having a continu- 
ous distribution of local slopes. In contrast, the ordered 
surfaces of sandpiles submitted to large tilt are globally 
smoother, despite being locally rougher, with a discrete 
distribution of local slopes. This preliminary indication 
of anomalous roughening Q in sandpile automata can be 
explained, tentatively, in terms of local slope variations. 

When disorder is present, the sandpile surface can 
'wander' over length scales corresponding to many 
columns by making small changes in local slope, thus 
adding to the surface width without compromising stabil- 
ity. However, in ordered sandpiles, the allowable changes 
in local slopes are comparatively large and the width can- 
not increase indefinitely without making the sandpile un- 
stable. Our observation of an abrupt decrease JjJ in the 
surface width of the sandpile at the dotted line separat- 
ing triangular and uphill avalanches in Fig. [l], supports 
this argument. 

Lastly, we provide a theoretical framework for the on- 
set of uphill avalanches in terms of directed percolation. 
Uphill avalanches occur when the bonds that transfer 
slope instabilities are connected in an upwards direction; 
it is worth noting that although the appropriate percola- 
tion threshold should be that of directed bond percolation 
(0.64), the pre-existence of downhill avalanches before 
this onset implies that the the backbend length |lC[ ] is 
infinite in the downward direction. Hence the per- 
colation threshold is reduced to 0.5, the undirected bond 
percolation threshold for square lattices. From Eq. 2, 
it can be shown that an instability at site (i,j) is trans- 
ferred uphill when the height decrease 5z(i,j) at that site 
satisfies one of the two following conditions : 

5z(i,j) >2(s c -s(i + l,j)) 
Sz(i,j)>2(s c -s(i,j+l)) (4) 

This can be used to write, in mean field, the compound 
probability T* for the uphill transfer of instability 

T T =p(5z = 2)p(s c -s< 1) 

+ p(6z = 1 + a)p(s c - s < -(1 + a)) 

+ p(5z — 2a)p(s c — s < a) (5) 

where the p's represent the probabilities of the events 
in parentheses. These probabilities can be computed di- 
rectly from our simulations, and we find that T T varies 
near-monotonically between 0.45 (for tilt angles where 
triangular avalanches are predominant) and 0.55 (for tilt 
angles where uphill avalanches are the norm) . We suggest 
that the transition to predominantly uphill avalanches 
occurs when = 0.5, i.e. when there is an infinite 
cluster of bonds connecting unstable sites uphill from the 



point of initiation. This is in rough agreement with our 
simulation results; better agreement can only be obtained 
by going beyond mean field and including fluctuations. 
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FIG. 1. A stability diagram for two-dimensional sandpiles 
with L = 32, a — 0.7, d- = 2 and d\ = 20. The measured 
mean slope < s > is plotted against the critical slope s c , which 
should be interpreted as an inverse tilt (see text). The crosses 
(x) represent steady-state stability, the full line represents the 
spontaneous flow threshold, and the broken line separates a 
region of triangular, predominantly downhill avalanches (A) 
from that of uphill avalanches (•) 



FIG. 2. Contour plots showing the frequency with which 
avalanches initiated near the centre of a sandpile cover its 
different regions. The piles have L = 50, a = 0.7, d- = 2, 
d\ = 20, s c = 2.05 and (a) s' c = 1.95 (b) s' c = 1.75. The 
contours correspond to 20, 40, 60, and 80 per cent coverage, 
and the images are built from ~ 3500 test events. 
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FIG. 3. Distribution functions for local slopes in disor- 
dered sandpiles with L = 50, a = 0.7, d- = 2, and d\ = 20. 
(a) A steady state pile with critical slope s c = 1.85. (b) A 
pile subjected to small tilt; s c — 2.05 and s' c = 1.85. (c) A 
pile subjected to large tilt; s c = 2.05 and s' c = 1.75. 



FIG. 4. The mean distance (•) between the centre of the 
avalanche and its point of initiation, and the average fraction 
of avalanche sites uphill from the initiation point (x) , plotted 
against the degree of tilt for a model sandpile with parameters 
L = 50, a = 0.7, d- = 2,d\= 20 and s c = 2.05. 
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